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Nanoparticles can influence the properties of polymer materials by a variety of mechanisms. With 
fullerene, carbon nanotube, and clay or graphene sheet nanocomposites in mind, we investigate how 
particle shape influences the melt shear viscosity rj and the tensile strength r, which we determine 
via molecular dynamics simulations. Our simulations of compact (icosahedral), tube or rod- like, 
and sheet-like model nanoparticles, all at a volume fraction ~ 0.05, indicate an order of magnitude 
increase in the viscosity n relative to the pure melt. This finding evidently can not be explained by 
continuum hydrodynamics and we provide evidence that the rj increase in our model nanocomposites 
has its origin in chain bridging between the nanoparticles. We find that this increase is the largest 
for the rod-like nanoparticles and least for the sheet-like nanoparticles. Curiously, the enhancements 
of rj and r exhibit opposite trends with increasing chain length iV and with particle shape anisotropy. 
O , Evidently, the concept of bridging chains alone cannot account for the increase in r and we suggest 

that the deformability or flexibility of the sheet nanoparticles contributes to nanocomposite strength 
and toughness by reducing the relative value of the Poisson ratio of the composite. We note the 
molecular dynamics simulations in the present work focus on the reference case where the modifi- 
cation of the melt structure associated with glass-formation and entanglement interactions should 
not be an issue. Since many applications require good particle dispersion, we also focus on the 
case where the polymer-particle interactions favor nanoparticle dispersion. Our simulations point 
to a substantial contribution of nanoparticle shape to both mechanical and processing properties of 
polymer nanocomposites. 
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I. INTRODUCTION 



Blends of polymers and nanoparticles, commonly called "polymer nanocomposites" (PNC), have garnered much 
attention due to the possibility of dramatic improvement of polymeric properties with the addition of a relatively 
J> ■ small fraction of nanoparticles [l|, 0, [H, 0, [E IE Q- Successfully making use of these materials depends upon a 
firm understanding of both their mechanical and flow properties. Numerous computational and theoretical studies 
have examined the clustering and network formation of nanoparticles and their effect on both the structural and 
Theological properties of PNCs @, |, EE El El El El EE EE El El [ll EH HI H, H H, HI] . The vast majority 
of these efforts have focused on nanoparticles that are either spherical, polyhedral or otherwise relatively symmetric, 
although there are some notable exceptio ns [l9lj20lj2ll[23l. In contrast, experiments have tended to emphasize highly 
asymmetric nanoparticles (27l. l2cl l29l. l3f]|. l3ll, I32L l33l. 134 l35l l36l. l37l. l38j . such as layered silicates or carbon nanotubes. 
It is generally appreciated that these highly asymmetric nanoparticles have the potential to be even more effective 
than spherical (or nearly spherical) nanoparticles in changing the properties of the polymer matrix to which they are 
added. In addition to the large enhancements in viscosity and shear modulus expected from continuum hydrodynamic 
and elasticity theories, extended nanoparticles can more easily form network structures both through direct interaction 
between the nanoparticles, or through chain bridging between the nanoparticles [IJ EE E3> > where a "bridging" 
chain is a chain in contact with at least two different nanoparticles. These non-continuum mechanisms are believed to 
, play a significant role in property enhancement, though the dominant mechanism depends on the properties considered, 
particle-polymer and particle-particle interactions, sample preparation, etc. 
. i-H ! Given that the majority of previous computational efforts have focused on symmetric nanoparticles, we wish to 
elucidate the role of nanoparticle shape in determining basic material properties, such as the viscosity 77, and material 
"strength", (i.e., breaking stress). Computer simulations are well suited to examine the role of nanoparticle shape, 
since it is possible to probe the effects of changing the shape without the alteration of any of the intermolecular 
interactions. In this way, the changes due to nanoparticle shape can be isolated from other effects. Such a task is 
complicated experimentally, since it is difficult to modify the shape of a nanoparticle without dramatically altering 
its intermolecular interactions. In this paper we evaluate the viscosity 77 and ultimate isotropic tensile strength r 
of model PNC systems with either (i) highly symmetric icosahedral nanoparticles (compact particles), (ii) elongated 
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rod-like nanoparticles, and (iii) sheet-like nanoparticles. These nanoparticles can be thought of as idealizations of 
common nanoparticles, such as gold nanoparticles and fullerenes (polyhedral), nanotubes and fibers, and nanoclay 
and graphene sheet materials, respectively. Our results are based on molecular dynamics (MD) computer simulations, 
using non-equilibrium methods to evaluate r\ (40l . l4lj , and exploiting the "inherent structure" formalism to determine 
t [42|, |43j]. We find that the rod- like nanoparticles give the largest enhancement to 77, which we correlate with the 
presence of chains that bridge between the nanoparticles. The sheet nanoparticles offer the weakest increase in 77, and 
correspondingly have the smallest fraction of bridging chains. For the ultimate isotropic strength r, we find opposite 
results: the sheets provide the greatest reinforcement, while the rods the least. For both of these properties, the 
property changes induced by the icosahedral nanoparticles fall between those of the extended nanoparticles. 

The present simulations are idealized mixtures of polymers and nanoparticles in which the polymer-nanoparticle 
interactions are highly favorable so as to promote nanoparticle dispersion. Moreover, we have chosen to work at 
relatively high temperature in order to avoid contributions to 77 from the complex physics of slowing dynamics 
approaching the glass transition. Previous work 0, [l(| has shown that polymer-surface interaction effects in this low 
temperature range can alter, and potentially dominate the nanocomposite properties. We also limit the range of chain 
length N studied to avoid effects of significant polymer entanglement. These limitations on interaction, temperature, 
and chain length are advantageous in order to develop a clear understanding of the origin of the observed changes 
in properties. Such a reference calculation provides a reference starting point to understand behavior when these 
constraints are relaxed. With this in mind, caution is needed when comparing these results with experimental data 
where these complicating additional factors may be present - along with other possible effects, such as crystallization 
or phase separation. 

We organize this paper as follows: in Section HH we describe the details of the model and method, focusing on 
the differences between the nanoparticle types used in each system. Section IIIII describes our investigation of the 
rheological properties of the nanocomposites, while Section IIVI considers the effects of shape on r. We conclude in 
Section El 



II. SIMULATION 



To directly compare to experiments, it is desirable to use as realistic a molecular model as possible. While a 
chemically accurate MD simulation is possible in principle, it is often more difficult to identify basic physical trends 
with such models. Such attempts at chemical realism are also demanding in terms of the computational times 
required, which restricts the class of problems which can be investigated. Coarse-grained models of polymeric materials 
provide a good compromise between the the opposing needs of realism and computational feasibility. Such models 
can reproduce qualitative experimental trends of nanocomposites, butprecise quantitative predictions cannot be 
expected [26[ . Building on nanocomposite models introduced before [1, [tj Qil 13 , we study a coarse-grained model 
for polymers [45|, [4|| and nanoparticles that allow us to consider PNC systems over a wide range of physically 
interesting systems. Here we consider several nanoparticle shapes built by connecting spherical force sites. 

We perform MD simulations of systems consisting of a small fraction of model nanoparticles in a dense polymer 
melt. For reference we also simulate the corresponding pure polymer melt. The polymers are modeled via the common 
"bead-spring" approach, where polymers are represented by chains of monomers (beads) connected by bond potentials 
(springs) [47[. All monomers interact via a modified Lennard Jones (LJ) potential 
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where is the distance between two monomers, where ep_p is the depth of the well of the LJ potential and a 
is the monomer size. The potential is truncated and shifted at r c — 2.5 a, so that the potential and force are 
continuous at the cutoff. Bonded monomers in a chain interact via a finitely extensible, non-linear elastic (FENE) 
spring potential [HI, l4t| 
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where k and Ro are adjustable parameters that have been chosen as in ref. [45l . |46|. Since we do not aim to study 
a specific polymer, we use reduced units where a = ep-p = m = 1 (m is the monomer mass). Length is defined in 
dimcnsionlcss units relative to er, time in units of a^Jraje and temperature T is expressed in units of e/kg, where ks 
is Boltzmann's constant. 
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We use three different types of nanoparticles for our calculations. Fig. [T] shows representative images of the 
nanoparticles. The first type of nanoparticle, an icosahedron, was previously studied in ref. [1, 0, [1(3] which focused 
on factors controlling nanoparticle dispersion and low temperature effects on transport for a similar polymer matrix, 
respectively. Practical realizations include fullerene particles, primary carbon black particles, quantum dots, and 
metal nanoparticle additives. The nanoparticle force sites interact with each other via an identical Vp_p given in 
Eq. (H]), with (7 = ep_p = m = 1. To maintain the icosahedral shape, the force site at each vertex is bonded to its 5 
nearest neighbors via a harmonic spring potential 



VWnj) = -nf M - lj (3) 

where k = 60 and r equals the minimum of the force-shifted LJ potential, approximately 2 1 / 6 . To further reinforce 
the icosahedral geometry, a central particle is bonded to the vertices with the same potential Eq. with the same 
value of n, but with a slightly smaller preferred bond length, r' = 1/4(10 + 2 v / 5) 1 ^ 2 f*o, the radius of the sphere 
circumscribed around an icosahedron. The resulting nanoparticles have some flexibility, but are largely rigid, thereby 
preserving their icosahedral shape. 

The second type of nanoparticle is a semiflexible rod represented by 10 LJ force sites with neighboring monomers 
bonded by the same Vfene as used for the polymers. We choose the shape to represent nanoparticles such as carbon 
nanotubes or nanofibers. Carbon nanotubes typically have lengths up to several /xm and a diameter of 1 nm to 
2 nm (single- walled tubes) or 2 nm to 25 nm (multiwalled tubes can have even larger diameters) [48| . Unfortunately, 
such large rods are not feasible to simulate, as the system size needed to avoid finite size effects exceeds current 
computational resources. As a compromise, we simulate rods with a length-to-diameter ratio of 10. There is an 
additional bond potential between nanoparticle force sites 

Min(0)=KHn(l + OOS0) (4) 

where 9 is the angle between three consecutive force sites. This imparts stiffness to the rod. We choose Ku a = 50, so 
that like the icosahedra, the rods have some flexibility, but are largely rigid. 

The third type of nanoparticle is square "sheet" comprised of 100 force sites. While the sheets are obviously 
different from the rods, the sheets have the same aspect ratio as the rods, where aspect ratio is defined by the 
ratio of the largest and smallest length scales. Aspect ratio is often considered to be one of the most important 
properties for anisotropic particles, apart from the interactions with the surrounding matrix. The sheets represent a 
coarse-grained model of clay-silicate nanoparticles (23|, and in the literature these objects are also termed "tethered 
membranes" 0, [13, 5H, H(| |4{|. Each monomer in the 10 x 10 array has non-bonded interactions described by the 
same LJ forces as the polymers. The monomers are also bonded to their 4 nearest neighbors via the same FENE 
bond potential of Eq. ([?]). The particles at edges and corners of the sheet are only bonded to 3 and 2 neighboring 
particles, respectively. Like the rods, the sheets are stiffened with a Vu n potential to prevent the nanoparticle from 
folding in on itself. The sheets also include a perpendicular bonding potential 

V^) = ^(0-^) 2 (5) 
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which limits distortions from a square geometry. Without the potential Vj_, the sheet can deform into a rhombus. As 
in ref. [23|, we choose Kii n = 10 and k± = 100 for the sheets. 

Thus far, we have not defined a nanoparticle-monomer potential. In previous work 1^. fic| . the same LJ potential of 
Eq. (TT]) was used for nanoparticle-monomer interactions, with the well depth ep_p replaced by 6n-p- The soft form 
of the attraction makes it relatively easy for monomers which are first and second nearest neighbors of a nanoparticle 
to exchange. As a result, the chains can "slide" relatively easily along the nanoparticle surface, thereby reducing 
any benefits of networks built by chain bridging 50]. To make monomer exchange at the surface less favorable, we 
use a "12-24" potential U12-24 rather than the standard "6-12" powers of Vl,t- Thus, the the nanoparticle-monomer 
interactions are of the form 
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where a — 1 and r c = 2.5 a. This potential has shorter range and stronger forces binding the first neighbor monomers 
to the nanoparticle. So that the total energy of the potential well, i.e. 4ir JJ C l/N-pdr, is the same as that of the LJ 
potential used in ref. [U[l(3], we choose cn-p = 3. In this way, the total potential energy is comparable to that used 
in ref. [jjl fiot with €n-p = 1.5, which led to well-dispersed nanoparticles in that study. Thus, we expect our systems 
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will be well-dispersed. However, we shall see this is not entirely the case for the systems with sheet nanoparticles. 
Figure [2] shows snapshots of a typical configuration for each system studied. 

We generate initial configurations using the same approach as in ref. [1(3], namely growing vacancies so that 
the nanoparticles can be accommodated. However, this process generates artificial initial configurations, and so we 
generate subsequent "seed" configurations by simulating at T — 5 where reorganization occurs on relatively short 
time scales. Depending on the nanoparticle type, we extract independent starting configurations every 10 5 to 10 6 
time steps. Subsequently, we make the necessary changes in density and chain length before cooling and relaxing at 
T = 2. A possible concern is that at T = 2, monomers will stick to the nanoparticle surface on a time scale that 
is long compared to the simulation, since £n-p = 3; however, we confirmed that monomers of chains exchange with 
the nanoparticle surface many times over the simulation. The equilibration time depends on the system type and 
is directly related to particle size. Pure systems relax relatively quickly, needing only 10 5 time steps. Icosahedral 
systems need roughly 10 6 time steps before reaching thermodynamic equilibrium. The rods and sheet nanoparticles 
are much larger, and thus diffuse more slowly. As a result, these systems require in excess of 5 x 10 6 time steps to 
reach equilibrium. Given that equilibration requires more than 10 6 time steps at T — 2, equilibration under conditions 
of entanglement or supercooling would be computationally prohibitive. 

We integrate the equations of motion via the reversible reference system propagator algorithm (rRESPA) , a multiple 
time step algorithm used to improve simulation speed We use a basic time step of 0.002 for a 3-cycle velocity 
Verlet version of rRESPA with forces divided into "fast" bonded (Vfene, Vharm, Vim, Vj_) and "slow" non-bonded 
(Vp_p,Vn-p) components. The temperature is controlled using the Nose-Hoover method where the adjustable "mass" 
of the thermostat is selected to match the intrinsic frequency obtained from theoretical calculations of a face-centered 
cubic LJ system [5lj |. 

To study rheological properties, we shear equilibrated configurations using the SSLOD equations of motion [5l| . 
integrated using the same rRESPA algorithm used for equilibrium simulation. The SLLOD method generates the 
velocity profile for Couette flow. If we choose the flow along the x-axis and the gradient of the flow along the y-axis, 
the shear rate dependent viscosity is given by 

V = ~^- (7) 



where (P X y) is the average of the x-y components of the pressure tensor (sometimes also called the stress tensor), and 
7 is the shear rate [4l[ . We limit our simulations to small enough shear rates to avoid potential instabilities associated 
with breaking FENE bonds in the simulation [52| . 

We choose a loading fraction <f> w 0.05, defined by the ratio of the number of nanoparticle force sites to the total 
number of system force sites; values of this order of magnitude are common in experimental studies. [30l l34l [35l . [36j 
Since all force sites have the same diameter, </> should also be roughly equal to the volume fraction. Due to the fact 
both polymers and nanoparticles consist of a discrete number of force sites, <j) varies slightly between systems. For the 
icosahedral systems = 0.0494, and for the rods and sheets <j> = 0.0505. The system size for all three nanoparticle 
systems is much larger than either the radius of gyration of the polymers, or the largest nanoparticle dimension, to 
avoid finite size effects. The number of chains Nc ranges from Nc = 10 in the smallest system up to Nq = 994 for 
the largest systems. For each system we examine chain lengths N — 10, 20, and 40. For this polymer model, the 
entanglement length » 30 to 40 [53|. Hence, only the longest chains studied may exhibit any effects of entanglement 
and the effects should be small in the present work. We reiterate that effects of dynamics due to proximity to the 
glass transition should not play a significant role at the relative high T of our systems. A summary of the system 
parameters can be found in Table [T] 



III. COMPOSITE RHEOLOGY 



A. Shear Viscosity r\ 

In this section, we evaluate the role of nanoparticle geometry on the magnitude of r\ and clarify the influence of chain 
bridging on r\ for a range of shear rates (7 = 0.005, 0.007, 0.01, and 0.02) and chain lengths. In the nanocomposite 
simulations of [8j] , 77 became nearly constant at 7 = 0.005 - in other words, we are near the Newtonian limit where 
77 approaches a stationary value for 7 — ► 0. Thus, we expect our lowest shear rates approach the Newtonian regime; 
in this sense, the shear rates we use are fairly low. However, if we consider 7 in physical units, these rates would be 
extremely high by experimental standards. Ideally, we would estimate r\ in the 7 = limit directly (for example, by 
using an Einstein or Green-Kubo relation), but the accurate evaluation of 77 using these methods is difficult [13T, l44j|. 

Figure [3] shows 77 for the three nanocomposites and pure polymer as a function of 7. These data show a clear 
change in 77 between the systems for all 7 simulated, with the rods having the largest 77, followed by the icosahedra, 
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sheets, and finally the bulk polymer. The decrease of 77 with increasing 7 is indicative of shear thinning. We note 
that 77 is not constant at the lowest 7, indicating we have not yet reached the Newtonian regime. Nonetheless, it is 
clear that the pure polymer has a viscosity approximately an order of magnitude smaller than any of the composites. 
Comparable differences in 77 have been observed in experiments (29l . [3l| . This demonstrates an enhancement of 77 
through the addition of nanoparticles, which is often a desired goal of adding nanoparticles to a polymer melt. While 
the figure only shows 77 for the N = 20 and p = 1.00, other chain lengths and densities indicate the same trend. 
Figured] (a) shows the effect of varying N on 77 at p = 1.00 and 7 = 0.007. An increase in 77 with increasing N is 
expected from basic polymer physics since the chain friction coefficient increases with N [47] . Interchain interactions 
and "entanglement" interactions enhance this rate of increase since the friction coefficient of each chain increases 
linearly with N. 

While it is clear that addition of nanoparticles increases 77 of the resulting composite, the reasons for this are not 
obvious. Since the pure polymer has a relatively low 77 it is reasonable to assume that the rigid nanoparticles themselves 
inherently raise 77 as in any suspension of particles in fluid matrix. To separate out this continuum hydrodynamic 
effect of the nanoparticles from the effects specifically due to polymer-nanoparticle interactions, we calculate intrinsic 
viscosity [77] for the individual nanoparticles. Specifically, continuum hydrodynamics provides an estimate of the 
incremental change of ?7 r through the addition of particles to a fluid, 

r ?r (^) = l + [r ? ]^ + 0(^) (8) 

where the velocity of the fluid is taken to be zero on the surface of the particle ("stick" boundary conditions). Such 
expansions are only quantitatively reliable for </> < 0.01 for nearly spherical particles; for larger <p, we expect eq. (|8]) 
to recover only qualitative variations in r] r . For reference, [77] = 5/2 for spheres in three dimensions [541 ]. It is also 
known that [77] is generally larger for isotropically oriented spherical shapes [55[ . Refs. 0, [U, [56[ have reviewed 
the modeling of polymer composites containing extended particles based on continuum mechanics. This work also 
considers the effect of particle clustering on 77 which is important under poor dispersion conditions. 

To compare the predictions of hydrodynamic theory to our simulation, we use the Zeno package to calculate [77] 
for our particles [571 ] . The computational method involves enclosing an arbitrary-shaped object within a sphere and 
launching random walks from the border or launch surface. The probing trajectories either hit the object or return to 
the launch surface. From these path integration calculations, [77] for the rigid particles can be estimated for objects of 
essentially arbitrary shape [58(. We determine [77] of the nanoparticles as an average over 100 different confirmations 
found in the equilibrium system as they adopt a range of shapes due to bond flexibility. 

Using Zeno, we find that the icosahedral nanoparticles have [77] = 3.93, the rods have [77] = 16.3, and the sheets 



have [77] = 10.8. The intrinsic viscosity of a regular icosahedron has been estimated to be 2.47 |58[. The discrepancy 
between this value and our icosahedral nanoparticles arises from the fact that our nanoparticles are not regular solids, 
but rather consist of a group of spherical force sites with icosahedral symmetry. We can use Eq. ([H]) to predict 77,- 
for the MD simulations, and find (for the N — 20 systems) r\ Y — 1.20 for the icosahedra, r\ r — 1.82 for the rods, 
and 77,. = 1.54 for the sheets. Such a prediction underestimates the N — 20 MD results by 92 % for the icosahedra, 
90 % for the rods, and 89 % for the sheets. Clearly this hydrodynamic estimate of ?7 r for our nanocomposites is not 
adequate, even accounting for the order-of-magnitude nature of the estimate based on eq. ©. Consistent with the 
full MD simulations, the rods have the largest value, but in contrast to the MD simulation, the order of the sheets 
and icosahedra are reversed. While the larger [77] value of the rods are generally consistent with a larger value of 77, 
the reversal of the sheet and icosahedra [77] values indicates that nanoparticle and/or polymer interactions play a role 
in the overall viscosity. Moreover, the continuum theory indicates that 77 r should be independent of N, while we see 
from Fig. [4] (b) that 77,- strongly depends on N. Hence we conclude that the continuum model of 77 does not provide 
an adequate explanation of the changes we see in our simulations. 

Evidently, we must also examine other contributions to the nanoparticle viscosity. One possibility that has recently 
been discussed (2l| is that slip rather than stick boundary conditions might provide a better des crip tion of nanoparticle 
boundary conditions. However, this would result in a decrease of the predicted value of ?7 r [22l.l44| relative to the stick 
case, which is in the wrong direction for explaining our results. At a molecular level, the boundary conditions are 
clearly neither stick nor slip, as particles have finite residence times at the surface, although it is difficult to determine 
hydrodynamic boundary conditions directly from molecular considerations. Nanoparticle-polymer interactions and 
nanoparticle clustering arc evidently factors that might cause the property enhancements we observe, and in the next 
section we focus on this possibility. 



B. Relationship of Viscosity to Fluid Structure 

We next turn to the structural results obtained from MD simulation to better quantify the relationship of polymer 
and nanoparticle structure to 77. We first calculate the fraction of nanoparticles in contact with other nanoparticles, 
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/n-n- Nanoparticles are said to be "in contact" if one or more comprising force sites of a nanoparticle are within the 
nearest neighbor distance of a force site belonging to another nanoparticle. The nearest-neighbor distance is defined 
by the first minimum in the radial particle distribution function g(f); the first minima is at r ~ 1.5 for all systems. 
Figure [5] shows /n-n as a function of TV for p = 1.00 and for both 7 = and 7 = 0.007. While it is not readily apparent 
from Fig. [2](c), the polymers "intercalate" between the sheets to some degree, and hence, the nanoparticles are not 
actually in direct contract with each other, resulting in a very low value of /n-n for the sheets. With the exception of 
the N — 40 equilibrium system, the sheets have the fewest nanoparticle-nanoparticle contacts. For the systems under 
steady shear, /n-n evidently stays fairly constant as N varies. While the order of the systems for /n-n is consistent 
with the trend for 77, the similarity between the rods and icosahedra suggests that the fraction of nanoparticles alone 
is not responsible for determining composite viscosity. Of note, is the large fraction of nanoparticle contacts for the 
N = 40 sheet system at equilibrium. More in-depth analysis shows that the sheets favor a clustered state at lower 
density, i.e. chains no longer fully intercalate between the sheets. Both the N = 10 and N = 20 sheet systems have 
intercalating polymers, while the N — 40 case does not. The tendency of the sheets to stack is due in part to entropic 
interactions that result from depletion attractions between the sheet; the depletion interactions arise due to physically 
adsorbed polymers on the sheet surface, much like the depletion interactions of polymer coated colloids. (59| This 
interaction is most pronounced for the sheets, presumably due to the large and relatively flat surface of a sheet, which 
more effectively reduces the entropy of chains at the surface. 

Given that nanoparticle-nanoparticle contacts do not appear to be the origin of the relative ordering, we next 
consider the role of polymer-nanoparticle interactions. To gauge the possible role of such interactions we calculate 
the fraction of polymer chains in direct contact with a nanoparticle, /n-p- Using the same nearest-neighbor criteria 
as used to determine nanoparticle-nanoparticle contacts, we compute /n-p as a function of N for p = 1.00 at 7 = 
and 7 = 0.007 in Fig. [5] We find that /n-p is approximately the same for both the quiescent and sheared systems. 
There are also clear trends: /n-p increases with increasing N and for every value of iV the system containing the rods 
have the largest value of /n-p, followed by the system containing the icosahedra, and finally the system containing 
the sheets. These same trends in 77 appear in Fig. |4] suggesting that the nanoparticle-polymer contacts indeed have a 
significant effect on relative magnitude of 77. 

The correlation of 77 with /n-p is consistent with the idea that bridging between nanoparticles play a major role in 
determining the rheological properties of polymer nanocomposites [ll|, [lg, [U H^|. Thus, we extend our analysis to 
test for a correlation between 77 and chain bridging. We define a "bridging chain" as a chain that is simultaneously 
in contact with two or more nanoparticles. Figure [7] shows /b as a function of N for p = 1.00. The trend in /b is 
consistent with both Jn-p and 77, showing an increase in Jb with increasing N as well as a clear ordering between 
systems for every chain length. Although not shown, the results seen in Fig. [7] occur for every value of p that we 
have simulated. Hence, the fraction of bridging chains seems to be a useful "order parameter" for characterizing 
the polymer-nanoparticle interactions. If the polymer-nanoparticle interactions were sufficiently small, the bridging 
would not be expected to play significant role, and it is likely that the continuum hydrodynamic approach would be 
applicable. 

Similar to the notion of bridging is the idea that the nanoparticles can act as transient cross-linkers that can lead 
to effects equivalent with the formation of higher molecular weight chains [111, [3!| . To test this idea, we define an 
"effective chain" as the collection of chains that are connected by nanoparticles. We then define an effective chain 
length iVoff as the mean mass of chains connected by the nanoparticles. Fig. [8] shows that N c g is almost an order of 
magnitude larger than N . Hence, we expect that the largest contribution to the 77 increase comes from this effect. 
Even in the simplest Rouse theory, an order of magnitude increase in N would lead to an order of magnitude increase 
in 7;. The formation of longer effective chains is expected to lead to entanglement interactions that would further 
amplify the 77 increase, as emphasize by ref. '6Cj. However, this entanglement contribution is hard to quantitatively 
interpret in the present context. 

IV. ISOTROPIC TENSILE STRENGTH r 
A. Calculation of r 

The ultimate isotropic tensile strength r of a material is defined as the maximum tension a homogeneously stretched 
material can sustain before fracture. While this definition can be directly pro bed experimentally, the situation is less 
straightforward and computationally accessible in an MD simulation. Ref. [431 ] has developed an approach to estimate 
r based on potential energy-landscape (PEL) formalism that is accessible by MD simulation. Although the PEL is 
complex and multidimensional, one can envision it as a series of energy minima connected by higher energy transition 
pathways. By definition, the minima, or inherent structures (IS) are mechanically stable, i.e. there is no net force 
at the minimum. The method of ref. (43| relates r to the maximum tension the IS can sustain, determined by an 



6 



appropriate mapping from equilibrium configurations on the PEL. In other words, r is the upper bound on the tension 
at the breaking point of an ideal glass state in the T — > limit. Here we evaluate r for the various possible nanoparticle 
geometries. As a cautionary note, we point out that this method has not been directly validated by comparison with 
experiments; thus conclusions drawn from these results should be considered tentative. Ref. [6l| discusses an alternate 
approach using the PEL to evaluate the elastic constants. 

To determine r, we must first generate the energy minimized configurations from the equilibrium configurations. 
For a system of N atoms in the canonical ensemble, the most commonly used mapping from equilibrium configurations 
to the IS takes each force site along the steepest decent path of the energy of the system. This configurational mapping 
procedure corresponds to the physical process of instantaneous cooling to the T — > limit to obtain an ideal glass 
with no kinetic energy. Therefore, by sampling thermally equilibrated configurations and minimizing their energies 
we can estimate the average inherent structure pressure, Pis, which will be negative when the system is under tension. 
The maximum tension then defines r. The value of Pis has been found to be weakly dependent on the T at which the 
sampling is performed. For example, while cyclopentane has a strongly temperature dependent structure, ref. [62j has 
shown that equilibrium T only weakly effects Pis- Since the structures in our systems do not display such dependence, 
we expect our results also to be independent of the starting equilibrium T. Thus, we can use the same systems at 
T = 2 that we used to determine r\ and still obtain reliable results for r - even though the starting configuration is a 
highly fluid state. This calculation of r should be an upper bound for the tensile strength that would be obtained for 
any finite T system, and hence is referred to the as the ultimate tensile strength. The energy minimization process 
eliminates the high frequency effects that would normally be associated with such high T states and thus the method 
refers to the strength of an ideal glass material having essentially no configurational entropy. 

Figure [H shows the Pis curves generated for the N = 20 variant of each system. We find that Pis decreases at 
large p, until a minimum value is reached. The density of the minimum, referred to as the "Sastry density" ps [42| . 
characterizes the density where the system fractures, and voids first start to appear in the minimized configurations. 
As p decreases below this point, Pis begins to rise. Since this is the greatest tension achievable, we have P^ AX = — r. 
For pure, icosahedral, and rod systems we find ps ~ 0.99, while for the sheets it is much lower ps ~ 0.92. Similarly, 
r is considerably larger for the sheets than for the other systems. We find that the increase in r for the sheets is 
not as sizable as observed experimentally 29] . This may be related to the fact that (experimentally) the chains often 
form stronger associations with layered silicates than in our simulations; additionally, the chain lengths we examine 
are small compared to experiments, and we will see that r increases with N. We note that at this chain length, the 
icosahedra and rods actually reduce the strength of the nanocomposite as compared to the pure melt. 

To quantify the chain length dependence of r, we plot the minimum of each Pis curve as a function of N in Fig. 1101 
For the pure polymers, we find that r decreases with increasing chain length. This is non-trivial, since one might 
naively expect longer chains to exhibit more interchain coupling. This is the same trend with chain length as observed 
in ref. [431 ] for n-alkanes. It is reassuring that simple bead-spring model leads to similar results but this does not help 
us build intuition about the physical meaning of these results. Reference (43j finds a maximum t for N = 3, and a 
tendency for pis to saturate near N — 8. Since the smallest chain we simulate is N = 10, we are above the regime 
where this feature occurs. Consistent with these facts, r decreases with N and ps is roughly independent of N. 

Most importantly, Fig. [TU] shows that the addition of the nanoparticles reverses the iV dependence of r when 
compared to the pure system. Thus, while the presence of icosahedral or rod-like nanoparticles decreases the material 
strength for most chain lengths studied, the trend of r is increasing, and if this continues, r will surpass the pure 
melt for all nanoparticle shapes at large enough N. Indeed, r for the icosahedra nanocomposite already exceeds that 
of the melt at N = 40. 



B. Relationship of Structure to Tensile Strength 

To better understand the predicted N dependence of r, we perform a parallel analysis to the structural analysis of 
shear runs discussed above. The analysis of Figs. [5][7] focused onp= 1.00. We relate the structure to r by examining 
the structure at ps- Evidently, ps differs for the sheets in relation to the other systems studied, and it is possible 
that p dependent changes in connectivity properties are responsible for the difference of r. Thus, we calculate the 
quantities /b, /n-n> and /n-p at ps for each system using the equilibrium and IS configurations at the Sastry density 
ps of the system. The results for the equilibrium and IS configurations show the same qualitative trends, so we will 
present only data for the IS. 

We first focus on /b because our investigations into the rheological properties suggested that bridging chains played 
a large role in determining rj. However, when we plot /b in Fig.[TT]for each system at ps, we see that the N dependence 
is the reverse of that seen for r, even though this quantity evidently follows 77. Even for the sheets which have a ps 
significantly lower value of than the other systems, the ordering of /b among the systems is the same as for /b in the 
case of the 77 calculations. The sheet composites have the smallest /b, followed by those with icosahedra, and those 
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with rods. While bridging chains do increase with increasing N, /b does not seem to be a major factor in the relative 
value of t between the composites. For example, the N — 10 nanocomposites with icosahedra and sheets have almost 
the same value of /b, yet they have significantly different values for r. 

The fewer bridging chains in the systems with square sheets suggests that the sheets may be clustered, and we find 
that this is indeed the case for N = 40 (Fig. [5] (a)). Thus, we plot /n-n in Fig. [T2l and find indeed that /n-n ~ 1 
for the sheets. The sheets prefer a stacked state at low p due to entropic interactions with polymers which leads to 
a reduced explicit energetic interaction with the surrounding polymers. In particular, the relative ordering of /n-n 
matches that of r. Naively, this would seem to suggest a potential correlation between /n-n and r. However, this 
apparent correlation is problematic. Firstly, the N dependence of r and /n-n ar e different, namely r monotonically 
increases, while /n-n is roughly constant. Secondly, if the nanoparticles interactions are the origin of the increased 
strength, then one might expect to find the largest effect when the nanoparticles' are well dispersed Q. This is not 
the case, since the clustered sheets give the largest effect. 

Ref. [43j demonstrated that the rupture of the system occurs in "weak spots." To confirm that the nanoparticles 
actually impart additional strength, we visually examined the location of fracture in our systems at ps- To find empty 
spaces that are at least the size of a particle, we discretize the system into a cubic lattice of overlapping spheres, each 
with a radius ry and a nearest neighbor separation d. We then identify spheres that do not contain any system force 
site. By adjusting the parameters ry and d we can ensure that the voids we find are at least large enough for a single 
particle. Since the force sites of both the nanoparticles and polymers have diameter 1, we must choose ry > 0.5 to 
have physical relevance. We find that values of ry = 0.6 and d — 0.3 work best for visualizing voids. 

Figure [T31 shows that the voids in the system occur in regions of pure polymer for the icosahedral system. The same 
is also true for the rods and square sheets. Quantitative analysis of the force sites within nearest-neighbor distance 
r = 1.5 of the voids confirmed this suggestion, with over 99 % of the resulting force sites belonging to polymer 
chains. Thus, while the nanoparticles impart an increased strength to the nanocomposite, the correlation to /n-n is 
not apparently causal. Evidently, there is something more subtle controlling the magnitude of r in the sheet filled 
nanocomposite that we have not yet identified. 

To complete the structural analysis we calculate /n-p for the IS. Figure [14] shows that the sheet nanoparticles are 
not in contact with many polymers, consistent with the fact that /n-n ~ 1- In fact, after minimization the difference 
in /n-p between the rods and the sheets is even larger than observed in the equilibrium configurations shown in Fig. [6] 
Figure [TJ] shows that the relative ordering of /n-p is opposite to that of r, hence the number of nanoparticle-polymer 
contacts alone also does not provide a good indicator of r. 

A potential clue to the increase of r of the sheet nanocomposite is the fact that ps is significantly smaller than 
for the other systems. This allows the sheet composites to undergo a greater deformation before fracture. This 
large increase in both the strength and toughness of the polymer matrix with incorporation of nanoparticles is 
reminiscent of the chan ges in the properties of natural and synthetic rubbers with the inclusion of carbon black and 
nanofiller additives [39l. l63l. [H, [65 1. Extraordinarily large increases in both the strength and toughness of materials 
have recently been observed in the case of exfoliated clay sheets dispersed in the polymer polyvinylidene fluoride 
(PVDF) [6Q|. Although some of this change is associated with the modification of the crystallization morphology 
by the clay nanoparticles, this does not explain in itself the observed toughening mechanism. It is known that 
nanofiller particles can behave as temporary cross-linking agents that impart viscoelastic characteristics to the fluid. 
Gersappe [39l ] further emphasizes the role of this transient network formation in impeding cavitation events that 
initiate material rupture and the potential significance of this phenomenon in understanding the nature of biological 
adhesives (abalone), fibers (spider silk), and shell material (nacre) [6?], HI]. Gersappe further suggests that the relative 
mobility of nanoparticles has a role on the toughening of materials. However, such reasoning is inconsistent with the 
fact that the sheet nanoparticles lead to the strongest material while being the least mobile of our additives. 

Recent work indicates that flexible sheet-like structures, such as studied here, are characterized by a negative Poisson 
ratio v (69l . [TOL f7l| , i. e. the material expands normal the the direction of stretching, as opposed to normal materials 
which expand in the same direction as the stretching. The formation of a composite material with negative Poisson 
ratio ( "auxetic" materials) is expected to lead to a reduction of the Poisson ratio of the composite as a whole [72|, [7|[ . 
A large reduction in v is expected to give rise to materials that are strong and fracture resistant (74|. Related to 
this, we observe in our simulations that the voids that form in the sheet-filled nanocomposites tend to be fewer and 
smaller than those for the other nanoparticles at ps- As we further reduce p for the sheet system, the voids grow 
more numerous as opposed to growing larger as in the other nanocomposites. Microvoid formation has also been 
established as a mechanism for toughening polymer materials [75| . Previous work devoted to understanding the high 
impact strength of polycarbonate and other glassy polymers has likewise emphasized the importance of large "free 
volume" within the polymer material as a necessary condition for large toughness [76l |77| . Based on our simulation 
results and previous observations, we tentatively suggest that the addition of these "springy" sheet materials reduces 
the effective Poisson ratio of our nanocomposites, and that the microvoid formation process that we observe is a 
manifestation of the non-uniformities in the elastic constants within the nanocomposite. In principle, the elastic 
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constants can be determined by the PEL approach [6l| , but this analysis will be deferred to a future work since it is 
rather involved. Our tentative interpretation of the predicted r variation in the sheet nanocomposites also suggests 
the need for better characterization of the geometric rigidity properties of the sheet nanoparticles, since these variables 
may be important for understating how such particles modify the stiffness and toughness of nanocomposites. 

V. CONCLUSION 

In this work we have focused on how nanoparticle shape influences the viscosity of polymer-nanoparticle melt 
mixtures at high temperature and the ultimate tensile strength (estimated from the PEL formalism) of polymer 
nanocomposites in the ideal glass state. Our results suggest that chain bridging between the nanoparticles can have a 
large effect on 77 of the mixture when the polymer-particle interactions are attractive, so that the nanoparticles disperse 
readily. In addition, there is a relatively weak increase in rj for sheet nanoparticle composites, which tend to cluster 
in our simulations, supporting the expectation that nanoparticle clustering diminishes the viscosity enhancement. 
(However, the formation of "open" or percolating fractal clusters may have the opposite effect on viscosity). Curiously, 
the tensile strength of the sheet nanocomposite is greatest, in spite of the sheet stacking. One of the most intriguing 
effects of the nanoparticles is that, regardless of shape, the dependence of the tensile strength on chain length for 
the nanocomposites is opposite to that for the pure polymer melt. We reiterate that our results for r rely on the 
PEL approach, and this approach should be carefully compared with experimental measurements to test the expected 
relationship between r from simulations and tensile strength found experimentally. 

The trends observed for the r are more difficult to understand than those for rj. There is evidently no clear-cut 
correlation of r with the formation of bridging chains or with the attractive particle-particle interactions. In the 
absence of such a correlation, we suggest that the nanoparticles effect the mechanical properties by modifying the 
ratio of the bulk and shear moduli in the low temperature nanocomposite state (i.e. the nanocomposite Poisson ratio). 
Recent work has noted that molecularly thin sheets have a negative Poisson ratio when they are flexible enough to 
crumple by thermal fluctuations. Such additives could reduce the Poison ratio of the material as a whole and provide 
a potential rationale for interpreting the increases in both the strength and toughness observed in our simulation of 
sheets dispersed in a polymer matrix, as well as in experiments on clay nanocomposites. 
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TABLE I: The table details all system variants simulated; listing chainlength N, number of chains Nc, loading fraction cj>, 
number of force sites per nanoparticle n, and number of nanoparticles nn- 
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FIG. 1: The three different types of nanoparticles used in simulation: (a) icosahedron, (b) rod, and (c) sheet. The nanoparticle 
force sites are rendered as spheres connected by cylinders representing FENE bonds. 




FIG. 2: Typical equilibrium configurations of each of the three nanocomposite systems: (a) icosahedra, (b) rods, and (c) sheets. 
In each image all nanoparticles are shown; for clarity of the figure, only 10 % of the polymers are shown. 
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FIG. 3: The viscosity r\ as a function of shear rate 7 for chain length TV = 20 at p = 1.00. The rods show the largest r), followed 
by the icosahedra, and lastly the square sheets. We show the statistical uncertainty of rj for each system at 7 = 0.005, where the 
fluctuations is largest, and hence represents an upper bound for the relative uncertainty for all rj calculations. The uncertainty 
is the result of "block averaging" [4j]. The inset includes the bulk polymer and shows that the pure system has much lower 
viscosity than any of the nanocomposites by roughly an order of magnitude. The lines in this and subsequent figures are drawn 
only as guide for the reader's eyes. 
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FIG. 4: The (a) viscosity r\ and (b) reduced viscosity r; r as a function of chain length N. All values are from systems at p = 1.00 
and 7 = 0.007. Chain length appears to have no effect on the ordering of rj amongst the systems. 



12 




FIG. 5: The fraction of nanoparticles in contact with other nanoparticles /n-n as a function of N for each of the composite 
systems at (a) 7 = and (b) 7 = 0.007. 
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FIG. 6: The fraction of polymer chains in direct contact with a nanoparticle /n-p shown as function of chain length TV for 
each of the three nanoparticle systems at (a) 7 = and (b) 7 = 0.007. 
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FIG. 7: The effect of chainlength N on the fraction of bridging chains /b in each of the three nanoparticle systems at p = 1.00 
for (a) 7 = and (b) 7 = 0.007. The order of /b with respect to nanoparticle type matches the trend in Fig. [3] for 77. 
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FIG. 8: The effective chainlength N e g as a function of -/V in each of the three composite systems at p = 1.00 for (a) 7 = and 
(b) 7 = 0.007. 




FIG. 9: The inherent structure pressure Pis as a function of p for the pure polymer and three nanocomposite systems at 
N = 20. The pure, icosahedron, and rod systems all have a Sastry density ps ~ 0.99, while for the sheets ps ~ 0.93. The 
uncertainty intervals ("error bars") represent the statistical uncertainty in our average for Pis from block averaging [4ll ]. 
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FIG. 10: Tensile strength r as a function of chain length N for the pure polymer and three nanocomposite systems. We find 
that as r increases as N increases. This is in contrast to the behavior of the pure polymer where r decreases with increasing N. 
Only for systems with the sheet nanoparticles — or the system with icosahedral nanoparticles at TV > 40 — does the addition 
of nanoparticles produce a net benefit relative to the bulk polymer. 
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FIG. 11: The effect of chainlength N on the fraction of bridging chains /b in each of the three composite systems calculated 
using the inherent structure configurations at p = ps. 
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FIG. 12: The fraction of nanoparticles in contact with other nanoparticles /n-n as a function of N for each of the composite 
systems calculated using the inherent structure configurations at p — ps. 




FIG. 13: Typical fractures for the IS of the icosahedral nanocomposite. The large contiguous blobs are the fractures. The 
surrounding spheres represent force sites in contact with the rupture. Note that nearly all sites in contact are polymers. 
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FIG. 14: The fraction of polymer chains in direct contact with a nanoparticle /n-p as function of chainlength N for each of 
the three nanoparticle systems calculated using the inherent structure configurations at p = ps. 
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